Cold and Hot Dark Matter from a Single Nonthermal Relic 
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The origin of dark matter in the universe may be scalar particles pro- 
duced by amplification of quantum fluctuations during a period of dilaton- 
driven inflation. We show, for the first time, that a single species of particles, 
depending on its mass and interactions, can be a source of both cold and 
hot dark matter simultaneously. Detection of such weakly interacting parti- 
cles with masses below a fraction of an eV presents a new challenge for dark 
matter searches. 



It is widely accepted that a substantial fraction of the total energy density in the universe 
is in the form of dark matter (DM). The composition and amount of DM is not known, but 
many hypothetical particles were proposed as candidates: WIMPs, axions, LSP's, massive 
neutrinos and more DM is classified into two types according to the velocities of 

particles at the beginning of the epoch of structure formation when the temperature of the 
universe was about 1 eV, cold dark matter (CDM) if the particles are non-relativistic and 
hot dark matter (HDM) if they are relativistic. The prevailing wisdom is that CDM and 
HDM originate from different species of particles, for example, axions which weigh a fraction 
of an eV for CDM, and a few eV neutrinos for HDM. DM particles are traditionally assumed 
to have a thermal distribution of velocities (see, however, ||). Since a thermal distribution 
is sharply peaked around a single velocity, if that velocity is relativistic, the amount of 
non-relativistic particles is extremely small and vice versa. 



We show, for the first time, that weakly interacting nonthermal relics produced by am- 
plification of quantum fluctuations during a period of dilaton-driven inflation can have an 
energy spectrum with two peaks, such that at the beginning of the epoch of structure for- 
mation some fraction of the particles are relativistic and some fraction are non-relativistic. 
The ratio and magnitude of energy densities in relativistic and nonrelativistic particles are 
determined by the mass of the particles and by the cosmological parameters of the model. 
The basic physics behind the appearance of a twin peak spectrum is the existence of two 
scales in the problem: the Planck scale, redshifted, and the mass. Our models suggest the 
possibility that DM in the universe might be composed of weakly interacting nonthermal 
relics and that it might be a mixture of CDM and HDM if the mass of the particles is lighter 
than a fraction of an eV. It is argued j| that CDM+HDM models could be just what is 
needed to explain the data on CMB anisotropies ||. Our models can provide the desired 
composition from a single source, as well as predict other observable consequences. 

We consider particle production in models of string cosmology which realize the pre- 
big-bang scenario [§,0. In this scenario the evolution of the universe starts from a state of 
very small curvature and coupling and then undergoes a long phase of dilaton-driven kinetic 
inflation and at some later time joins smoothly standard radiation dominated cosmological 
evolution, thus giving rise to a singularity free inflationary cosmology. Particles are produced 
during the period of dilaton-driven inflation by the standard mechanism of amplification of 
quantum fluctuations [§]. 

In the simplified model of background evolution we adopt, the evolution of the universe 
is divided into four distinct phases with specific (conformal) time dependence of the scale 
factor of the universe a(r]) and the dilaton 0(77). The first phase is a long dilaton-driven 
inflationary phase, the second phase is a high curvature string phase of otherwise unknown 
properties, followed by ordinary Friedman- Robert son- Walker (FRW) radiation dominated 
(RD) evolution and then a standard FRW matter dominated (MD) evolution. We assume 
throughout an isotropic and homogeneous four dimensional flat universe, described by a 
FRW metric. All other scalar fields are assumed to have a trivial vacuum expectation value 
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during the inflationary phase. 

During the dilaton-driven inflationary phase rj < rj s , both scale factor and coupling 
are growing as powers a{rf) = a s anc l e Hv) = e 4>s (^-\, where a and (3 are negative. 

The dilaton-driven phase is expected to last until curvatures reach the string scale and 
the background solution starts to deviate substantially from the lowest order solution. For 
ideas about how this may come about see ||. The string phase lasts while r) s < 77 < 771. 
We assume that curvature stays high during the string phase. As in we assume that 



the string phase ends when curvature reaches the string scale M s . We parametrize our 
ignorance about the string phase background, as in [[□]] , by the ratios of the scale factor and 
the string coupling g(rj) = e^* 7 -*/ 2 , at the beginning and end of the string phase z$ = clx/cls 
and gi/gsy where g\ = e^ Vl ^ 2 and g$ = e^ Vs ^ 2 , where as = a(r) s ) and (fis = (j>(Vs)- We take 
the parameters to be in a range we consider reasonable. For example, z$ could be in the 
range 1 < zs < e 45 ~ 10 20 , to allow a large part of the observed universe to originate in 
the dilaton-driven phase, and gi/gs > 1, assuming that the coupling continues to increase 
during the string phase and 1CT 3 ^ g\ ^ 1CT 1 to agree with the expected range of string mass 
(see e.g. []IU|1). Some other useful quantities that we will need are U\, the frequency today, 
corresponding to the end of the string phase, estimated in TOJ to be uj\ ~ 10 10 Hz, and the 



frequency us = u\/zs, the frequency today corresponding to the end of the dilaton-driven 
phase. In the RD phase and MD phase are assumed to follow the string phase, the dilaton 
is taken to be strictly constant, frozen at its value today. 

We have computed the spectrum of produced particles for the models described previ- 
ously fl2f . We have solved a linear perturbation equation, x'i + + M 2 a 2 — ^ Xk = 
where s(r)) = a(r]) m e l ^ v ^ 2 = a™e^ s//2 (jf-^ ^ " , imposing initial conditions corresponding 
to normalized vacuum fluctuations. The parameter m depends on the spin of the particle 
and I depends on its coupling to the dilaton. The perturbation first "exit the horizon" when 
krj ~ 1, when curvatures become larger than their wavelength, then they are "frozen" out- 
side the horizon when krj < 1, and then "reenter the horizon" at r] re . A duality symmetry 
|T3| exchanging the perturbation and its conjugate momentum can be used to follow the 



evolution of the perturbation equations for times at which the background evolution is not 
known precisely. To read off the spectrum for a more general case of background evolution 
that we consider here, the only required substitution in the results of |12[ is n s = \ — am — ^l. 



see also ||14|). Similar calculations have been performed by several groups and the results 



agree [|15| HT?||, whenever a comparison was possible. 

We will consider weakly interacting scalar particles, abundant in string theory and su- 
pergravity. For scalar fields m = 1, and we will consider for concreteness the following values 
for /, I = —1,0,1 corresponding, respectively, to moduli (including the dilaton), Ramond- 
Ramond axions, and Neveu-Schwartz axions. We will assume that the produced particles 
interact so weakly, that their interactions and decay are not sufficient to alter the primordial 
spectrum substantially. The particles we have in mind have typically gravitational strength 
interactions, which is definitely weak enough to satisfy our assumption, and masses below a 
fraction of an eV. 

A typical spectrum of a light scalar may be divided, at a given time, into three frequency 
regions: i) The massless region, ujs > uj > M. In this region particles are relativistic. ii) The 
"false" massive region, M > uj > uj m , where uj m = iO\ {M/Ms) 1 / 2 . In this region particles are 
nonrelativistic, but have reentered the horizon as relativistic modes. Hi) The "real" massive 
region uj m > uj. In this region particles are non-relativistic, and have reentered the horizon 
as non-relativistic modes. Note that physical frequencies redshift as the universe expands, 
and therefore boundaries of regions change in time. 

Different spectral shapes may result depending on parameters. For the interesting case, 
the spectrum increases with uj in the massless region, decreases in the "false" massive region 
and increases in the "real" massive region. In this case, the energy density in relativistic 
particles Qrel — jfp- = and that of nonrelativistic particles is given by Qnr — 
(uj = uJ m ). For this case the spectrum takes the following approximate form at the 
beginning of the epoch of structure formation rj eq , 
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where x = 2 + 2a + and A/" is a numerical factor, estimated in [|l^], which we will set to 
unity in what follows. 

We have ignored, so far, particles produced during the string phase, since that phase is at 
the moment less well understood. If the spectrum decreases there, then our approximations 
remain valid. If, however, the spectrum increases there, then a good approximate relation 
to use would be VL REL ~ (lo = uj{) . Otherwise more parameters describing the string 
phase should be added. Since this is not relevant to our main point, we choose not to do so. 

We impose constraints on the spectrum, and show that it is possible to satisfy all of them 
by giving a specific example. First we require that the desired spectral shape is obtained 

0<ac<l. (2) 

Then we require that some relativistic DM particles and some nonrelativistic DM particles 
exist at r] eq , 

M < UJ S (Veq), U m > U eq , (3) 

where u eq = ffl2 ^ ) ^j(Veg)- We also require that the energy density in produced particles is 
the main source of energy density in the universe and that approximately equal amounts 
exists (later we will see that this condition can be significantly relaxed). 

^massive i^eq) — 1) ^massless ij^eq) — 1- (4) 

Assuming that the level of fluctuations at very large scales is well below the CMB level (as 
we show later this must be the case), we require 

-fL ( W = uj ) < 10- 5 . (5) 
a In a; 

A typical interesting spectrum is shown in Figure [I], 
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FIG. 1. A typical twin peaked spectrum. The dashed lines represents two possible string phase 
spectra. 



The constraints on the spectral parameters can be translated into constraints on cosmo- 
logical parameters. The range of powers a and (3 allowed by condition (||D is summarized in 
Table |, 

TABLE I. Allowed range of cosmologies 



I = -1 



I = 



I = 1 
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-1 < a < 
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Condition (f?) leads to 



M>MA — 



Condition M) leads to 
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(6) 



(7) 



and 

Since cu m = ui(M/M s ) 1 ^ 2 and cog = coiZg 1 , we obtain by substituting (|7|) into (^|) 
(€) (1+ " )/2 - L Condition © leads to g\ g (g)" (ff ) 1/2 < 10~ 5 , which, using 

( \ 1/2 / \ x 

~tn) (wi) z s < 10~ 5 - We therefore obtain the following two conditions, 



M > 10 10 /*M< ^ 



o>i/ 



which, using the following numerical values for parameters P,|10|, uj ~ h x 10 30 eV, o;i 



10 _1 eV,Af fl =0i x 10 28 eV (where 0.5 < h < 0.8 and 10~ 3 < #i < 10" 1 ), become 

M > h 2 gi x 10 10/a; - 30 eV, 

2 S < / l -l/*-^-l/* 10 29-8/«-5/^_ (1Q) 

We now present a concrete example which serves to demonstrate that a reasonable range 
of parameters exists, in which all conditions are naturally satisfied. We look at axions 
(7 = 1) with masses below .1 eV for a cosmological model described in 0. In this model 
a = —2/(d + n + 3) = —1/6 and ft = —Ad/ (d + n + 3) = — 1, for d = 3, n = 6, which 
is in the range specified in Table |[ For this specific model x = 2/3, ^rbl — g\ {^) > 

"™^te) 2 £te)~ 1/3 - 

The ratio of relativistic to nonrelativistic energy densities is given by = 
(ff 2 -) , and since o; m = c^i (M/M s ) 1//2 we obtain 

Q>REL _ Uh ^-2/3 fM s ^ 5/6 



M 



From conditions (|10D we obtain conditions on M and zs~. 



M > h 2 gi x 10" 15 eV, 
1< z s < 6 x 10 5 h~ 5/2 g± 3/2 (12) 



Taking 10~ 10 eV < M < 10~ 2 eV, for which condition (||) is comfortably satisfied, we observe 
that the ratio Qrel '■ &>nr can vary in a range from well above unity to well below unity. 
For example, choosing g\ = lO^ 1 and g$ = 10~ 2 if the axion's mass is 10 _10 eV, and for 
zs~2x 10 4 we get VLrel '■ &nr =1:1 with both energy densities being near critical, and 
if we choose g\ = 10 _1 and g$ = 3 x 10 -2 , making Qrel — 1> and if z$ ~ 10 6 we obtain 
VLrel '■ &nr = 1 : 10, with Qnr — 1, and just for fun, the preferred ratio of 70% CDM 
and 20% HDM f|] is obtained for the same values of couplings, if Zs ~ 2 x 10 5 . Note that 
the ratio fly} ) depends on z s 

2/3 M -5/ 6; go the 

previous examples correspond to a range of 

allowed values. 

We now show that it is not possible to obtain HDM+CDM and produce the CMB 
fluctuations by the same field, accepting all conditions on the spectrum. If the spectral 
amplitude is indeed at the level of the observed CMB fluctuations of about 10~ 5 , then the 
spectral slope is also constrained by the data ||. By parametrizing the slope as dfl/d In uj ~ 
deduce that n can be identified with the tilt parameter of the spectrum fllS] 



which is constrained by the data to be in the range 0.8 < n < 1.4, which corresponds to 
—0.1 < x < 0.2. We have required a positive slope, so the available range is just < x < 0.2. 
Now it is a simple exercise to calculate the energy density at u m and see that it is much 
below unity, dO / d In uj{uj m )/ dVt/ d In uj{uj ) = (^)* Since ^ = ^(^) V2 < 10 16 , and 
since x < 0.2, we can estimate that this ratio cannot naturally be above 10 3 , making the 
energy density at uj m much less than unity in most of parameter space. Perhaps, by tuning 
and forcing parameters into corners, it is possible to find an artificial example, but we will 
not pursue this possibility. 

We further show that, in the range of parameters we are interested in, it is not possible 
that one scalar field provides the required CMB fluctuations and a different scalar field 
provides HDM+CDM in a single cosmology. Since the two fields must be different If ^ Ihc, 
where the subscript / denotes a field that is supposed to produce the CMB anisotropy and the 
subscript HC denotes the field that is supposed to produce HDM+CDM. We require ^Irel — 
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2 (^)' 2lHC * 1, and Q f REL ~ ^ (f ) 2h < 1. Therefore, ^ = (ja ) a ' > 1, so for 



^ < 1 we have to have If < Ihc and from Xf = 2 + 2a + lff3 and xhc = 2 + 2a + lucfi (recall 
that a < and /3 < ) we obtain Xf > xhc- The conclusion is that the only possibility 
(for values of I, —1, 0, 1) is < xhc < %f < 0.2. But then Q HC dominates also at the lowest 
frequency, in contradiction to our assumptions. Therefore, either CMB fluctuations are of 



different origin [|T^,|T^], or DM source is different. 

If DM in the universe is indeed made of light particles with gravitational strength in- 
teractions its detection in current direct searches is extremely difficult, and will probably 
require new methods and ideas. 

Finally, we note that the models we have described have predictions and consequences 
other than DM. Additional particles such as gravitons [11] get produced and provide addi- 



tional possible experimental and observational signatures. 
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